Regime interpretation of anomalous vortex dynamics in 2D superconductors 
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Low -frequency dynamic impedance {a^^{Lj,T) = (cri + i(T2)~^) measurements on Josephson junc- 
tion arrays with finite vortex screening length ^, found that cri ~ |logaj|, ~ constant. This 
imphes anomalously sluggish vortex mobilities ^v{oj) ^ crj"^, and is in conflict with general dynam- 
ical scaling expressions that yield, for low-tj, ai and a2 0. We calculate : a) a{u!,T) by 
real-space vortex scaling; b) /iy(tj) using Mori's formalism for a screened Coulomb gas. We find, 
in addition to the usual critical (large-tj) and hydrodynamic (low-oj) regimes, a new intermediate- 
frequency scaling regime into which the experimental data fall. This resolves the above mentioned 
confiict and makes explicit predictions for the scaling form of a{ui,T), testable in SNS and SIS 
arrays. 

PACS numbers: 74.50.4-r, 05.90.+m, 74.60. Ge 



The dynamic conductivity of superconductors 
a{uj,T) = (Ji + ia2 = lale**^", including high-Tc mate- 
rials and Josephson junction arrays (JJA), has been the 
focus of much recent interest . Dynamical scaling 

forms, a = ^^+^~'^S'±(F~^), (1)^ — ^{Y~^) were proposed 
by Fisher et al. ||] and Dorsey , where z is the dynamic 
exponent, and Y ~ The results apply also 

for the vortex- unbinding Kosterlitz-Thouless (KT) tran- 
sition in d = 2 ^ being the vortex screening length. 
The scaling functions and 5*+ {S-) have well-defined 
limits in the hydrodynamic (F ^ 1) and critical {Y ^ 1) 
regimes, e.g. ci —^ CT2 — » in the Y —>■ oo, dc limit. 
2D JJA's 1^ are clean, controllable 2D superconductors, 
and should be ideal systems to display the universal dy- 
namic scaling behaviour and limits. 
Remarkably, however, low-frequency dynamic impedance 
(cr~^) measurements on 2D SNS triangular-lattice 
JJA's with a (field-induced) vortex screening length ^, 
find fJi ^ I lna;|, (72 ^ constant, in conflict with dynamic 
scaling limits. This also implies anomalously sluggish 
low-frequency vortex mobilities, ^ erf ^ ~ 1/1 In'^l 

going to zero for w — > 0. cr is related to the dynamic 
dielectric function e(w) : a ^ lY/eiuj). Surprisingly, 
Minnhagen's phenomenology (MP) for e{Lo), described 
below, and related simulations support this anomalous 
behaviour, but understanding the apparent breakdown 
of dynamical scaling, in the very arena where one might 
expect its clear verification, is of central importance. 
In this Letter, we reconcile these results, by a "regime 
interpretation" |^ defined by the ratio Y — {r^^/^Y of 
the (squares of the) frequency-dependent diffusive probe 
length 1^ = ^JT\Ju} and the screening length ^. Here, 
Fq is a junction-determined phase diffusion rate, the 
lattice constant is unity, and we consider weak screen- 
ing and probes over several lattice constants : ^ ^ 1, 
» 1. We a) recalculate a{uj,T) by a real-space scal- 
ing 0, with an improved treatment of intermediate- 



scale screening; b) evaluate the vortex mobility /iy us- 
ing Mori's formalism for a screened Coulomb gas 
Three probe-scale regimes emerge. I) Probing free- 
vortex scales ("low" frequencies) ^ ^, in a "hy- 
drodynamic" region, Drude behaviour with the correct 
dc conductivity, ai{uj 0,r) ^ is recovered. 
II) At intermediate scales/frequencies, < ^, in a new 
" precritical" region, MP-like behaviour, cri ~ |loga;| is 
found. HI) Probing bound pair-scales ("high" frequen- 
cies) <C ^ in a "critical" region extending from just 
above Tkt to T = 0, a scale-dependent vortex damp- 
ing ^ (Ti ^ (r(^/lncj)^ ^ (a;(lna;)^)~^ is found, corre- 
sponding to large pairs moving in a logarithmically in- 
teracting viscous medium of smaller pairs. The results 
of a; ^ 0, T ^ T^j^, thus depend on the order of the 
limits. The ratio = ai(w,T)/a2(uj,T) = cot 4>a- at 
Y — 1 interpolates between Drude {Ra = 1) and MP 
(i?cr = 2/7r) signatures as w increases from zero, or 
T increases from T^j,. As a satisfying byproduct of cal- 
culation a), the MP-like expressions emerge as approxi- 
mations to the cr/^^ scaling function, valid in regime II. 
Calculation b) demonstrates that the general results are 
independent of the details of JJA dynamics, and depend 
only on Coulomb gas screening properties. Local spin- 
wave damping mechanisms specific to SNS arrays [pd], 12 
could play an additional role in producing anomalous be- 
haviour, widening regime II. But both SNS and SIS ar- 
rays should show all three regimes in principle, with dif- 
ferent relative sizes of regimes I and II, coming from very 
general considerations. 

Two different physical circumstances yield a non-zero 
free vortex density, and thus finite ^ : i) For zero external 



flux, and T > Tkt, = ^+\T) - 



^0 



above transition, T > Tkt] while = for T < Tkt- 
ii) Flux-induced vortices of concentration / ^ 1 Q (too 
dilute to form a stable lattice) can form a one-component 
plasma with a screening length ^, given by the Debye ex- 
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pression ^"^ = = (47r2//T)i/2 ^ 0, correspond- 

ing to "above transition" for any T. 
We now sketch the MP ideas 0, originally devel- 
oped to describe a{uj,T) structures at T = T^^ > 
Tkt, where ^+(T„) = r^^. The zero-wavevector con- 
ductivity a{uj,T) is proportional to the correspond- 
ing (inverse) dielectric constant: (t{uj,T)/(ToKo£^'^ — 
iY [ev{k = 0,uj,T)] , where co is a conductivity scale 
and Kq the bare vortex coupling. In MP, the real 
part 3?(£y(fc = 0,w, r)~^) of this zero-wave-vector dy- 
namic function is approximated by the zero-frequency 
static function, ^ (^ev{k,uj — 0,T)~^^ , evaluated at 
the probe scale, k — rj^. The imaginary part, 
S {sv{k = 0, w, is found from the Kramers-Kronig 

(KK) relations, that produce a InF dependence. Thus 
0, with sv{k,uj = 0, r)-i = fcV(fc2 + ^-2) : 



CT2 



Y 
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aoK„e Y + 1' 



2 Y^ InY 



aoKoe ttY^-1 



(1) 



At y = 1, Ra = 2/7r (i.e. (per — arctan (7r/2)), an 
MP signature. The dynamical scaling limits in the 
Y ^ I critical regime, both above and below Tkt, are 
ai (T2 ^ l/oj (independent of ^), $0- = 7r/2. In the 
y > 1 hydrodynamic regime, for T > Tkt [T < Tkt), 
one finds cti ~ 0-2 '^-^ 0, <i>cr '--^ (cri ^ S{uj), (72 ^ l/i^)- 
Eqn. (1) for MP has very different limits, however : 
CTi - £,'^Y^\\nYl (72 - l/tjforF < 1; and cri ~ C^llnFl, 
02 ^ i'^ iov Y :^ I. (Note that the F > 1 hmit, for 
^ = ■^+(r) fixed, ti> — > 0, implies infinite dc conductiv- 
ity, above Tkt)- We now outline our two complemen- 
tary calculations, with details elsewhere showing that 
MP-like behaviour emerges in an intermediate regime II, 
rather than in the scaling form regimes I, III. 



(A) REAL-SPACE VORTEX SCALING AND 

(j{io,T) 

The total (dimensionless) JJA bond current /*°* {f^—x, 
y directions) is a sum of Josephson or super- (~ 
sin A^^i), phase-slip or normal- (~ A^^^) and noise- cur- 
rents, and is conserved at every 2D lattice site i. If we 
ignore capacitive charge build-up on grains, 

iy^\A,A - A^S)) + U^it)] = 0. (2) 

Here, i^o = {2eRjIj/h)f = FqT, f-^ = {hIj/2ekBT), 
and /j, Rj are the junction critical current and RSJ 
model effective shunt resistance pO| , p^ , for the SNS/SIS 
array. The random noise current obeys {ffir{t)ffj.'r'{t')) — 
{2/h'o)5^fj_i5rr'S{t — t'). The JJA grain phases are — tt < 
9i < TT, and the external transverse vector potential 
A^i{t) = A^ii(u;)e~*"* is weak. Inverting the Laplacian 



A^ on 0i, the Langevin dynamics equation for the phase 
is MM : 



djm 



Fr>{t) 



(3) 



where f3H = -~ J2t,,r cos (A^^^ - Af,{t)), €„' = 
Grr is the 2D lattice Green's function (with singular part 
subtracted), and (^Fr{t)Fr' (t')^ = 2voGrr'S{t - t'). 
The dynamic conductivity calculation yields 
a^r.fj.'r' = o'p^r.fi'r' + S'^r.fj.'r'- Here, (T is the usual su- 
perfluid response, that at long wavelengths is (ct/ctq) = 
TTKoo^oS{^)+iKoa(ro/uj). With Op denoting an average 
with weight Pq — e^''^'^^^~"^(*)), a can be written as : 



T- 



sinA„0re~^''*sinA„/6'r' ) dt. 



(4) 



As before |lC|], we extract vortices by a dual transform, 
do a gaussian truncation on spin waves, and find that 
the effect of the Fokker-Planck "propagator" e~^°* is to 
produce a correlation angle decay as e"'""*. The correla- 
tion (T^r./^'r ' can be expressed as the projection (through 
derivatives) of the vortex partition " generating function" , 
with separated "test charges" at and /i'r'. Doing the 
time integral in Eqn. (4), the e~^°* factors lead to a "test- 
charge" (ia;/ro)/(l — iti-'/Fo) at /i'r' in the partition func- 
tion. The logarithmic potential, A'^Uo{R/a) — +2Tr6j^f^ 
and dipolar (^^^ — 0) scaling equations Q can be gener- 
alized 1^ to include weak (^~^ <^ 1) monopolar screening 
of a, a + da dipole-binding, approximated by a potential 
A^U{R/a) = +2iTgi8j^Q. Real-space integration of pairs 
of separation a, a -I- da can then be done, as usual pro- 
ducing the renormalized coupling Ki that obeys KT scal- 
ing equations. Since vortex damping is across the junc- 
tions in the array, it is scale-dependent, Fq — > F; = Fg/a^ 
[z — 2), where a = e', so the frequency-dependent test 
charge is {iujo^ /Tq) / {1 — iuia^ /Tq) After projection, 
this provides a dynamic Drude factor, at scale a = e', 
that weights the incremental, (static) scaling contribu- 
tions, d{Kigi). The KK relations are thus automatically 
satisfied. 

With a partial integration, the long-wavelength conduc- 
tivity aito) is then, finally, an integral over all pair con- 
tributions, with a range of length (and time) scales : 



= Y 



dlKigi 



ia/r^)^ 



dl 1 — i{a/r^Y^ 



(5) 



For a <C 1, 5/ — 1, and for a » 1, is the Debye 
dielectric constant / {q^ -f ^^) at a scale q ^ : gi ~ 
(1 -I- (a/^)^)^^. The dominant monopole effect on cr{Ljj) 
is the explicit gi factor, representing a scale-dependent 
reduction of far-off fields, as seen by a, a -\- da dipoles. 
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There is a smooth cross-over from dipolar (a <C ^, regime 
III) to monopolar (a ^ ^, regime I) screening with probes 
a, and with mutual (a ^, regime II) screening 
in between. Previously, we had matched regime I/III 
behaviour directly |10|, effectively taking gi to be a step 
function, and suppressing the intermediate regime. For 
external free-vortex screening — ^_d(/)), gi = 1/(1 + 
(o/C)^) throughout. 

Changing variables in Eqn. (5), a/i^ — > a we see a{uj, T) is 
a function of Y, with only logarithmic deviations ~ = 
In^, from the limits of the integral and in Ki Ki^i^. 
For ^ cx £,D{f) ^ v7i this implies a quasi-universality 
in F ~ //w, as found The imaginary part 02 of 
Eqn. (5) has a function peaked at a^/f^ = F in square 
brackets, multiplying a roll-off function. By rapid roll-off 
and sharp-peaking estimates, the total 02 is estimated as 
{l^ = Inr^) : 



C2 



Y 



Y- 



1 + Y- 



Y- 



l + Y 



(6) 



in the regimes I, II, III respectively, 01 Y ^ 1, Y < 1, 
y <C 1. Note in regime III, the factor cancels, and 
pof CT2/cro ~ Ki^/uj with the correct superfluid kinetic 
inductance limit Koc/i-^, for T < Tkt, uj 0. The real 
part of the total conductivity, apart from the (5(w) term, 
(using KK relations in regime II, where the integral is 
harder to estimate) is : 



Ki^ 2Ki^ InY 



1 dKi^ arctan Y 



i + r-2'7ri-r- 



2 dL 



Y- 



in regimes I, II, III respectively. The cr/^^ results of 
Eqns. (6) and (7) agree with the scahng limits ||J|] of 
S±(Y~^) and (j)cr{Y) in regimes I, III. In regime II, with 
dipolar screening neglected {Ki^^ Kq), oj^ is of the 
MP form, Eqn. (1), . With g if/ - if 00 + in regime 
III, c7i/cro ~ (ro/w)/Z2 for aU T < Tkt reflecting _the 
KT "critical line", in the dynamics; the phase angle 
4>cr = arctan (a2/cri) ^ 7r/2, as w ^ 0. 



(B) COULOMB GAS VORTEX DYNAMICS 

It is important to directly calculate the Coulomb gas 
vortex "charge" mobility ||ll| (for ■^^^ 7^ 0) in a way that 
is manifestly independent of the details ||l^ of JJA dy- 
namics, but shows the two frequency regimes I and II 
(since dipolar screening is not included, regime III will 
not appear). The overdamped equation of motion for a 

charge -f 1 is TyRi = — ^j^i ejVV{Ri — Rj), where the 
potential between charges (ej = ±1) in Fourier space is 
^io) = {q'^)^^ ^-nd Py is a friction coefficient. We use 
Mori's formalism |0,n3|, to relate /iy(w) to the corre- 



lation function $pp((f, w) for the vortex charge density 



P(^) = Si eiS{R— fi). The inverse mobility of a given 
particle or effective viscosity function is the sum of the 
bare friction coefficient and a contribution that is related 
to the forces from all other particles : 



) = Ty 


l + (fcBT)-l^|gy(<7)|^<i>pp(q,^) . 
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(8) 
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FIG. 1. Real and imaginary parts of the conductivity 
o"/o"oC^ versus y = Fo/o'C^, in a linear-log plot, from Eqn. (5). 
Here ^ is the vortex screening length and (Fo/a;)^''^ the probe 
length r^. The x symbols are experimental data points [4]. 
Inset ; regimes I, II, III in a schematic frequency-temperature 
diagram, with = C+^(r) defining the II/I boundary. 



(7) In order to evaluate <I'pp((f, oj) we make the usual ap- 



proximations 1^ (neglect of the "Mori projector" and 
decoupling higher correlations in terms of p(ci) and 71(17), 
the number density Fourier component). One obtains : 



(9) 



Here, Sp{q) is the static (charge) structure factor. Eqns. 
(8) and (9) determine ^v{^) self-consistently, but we 
solve to leading order, replacing /iy^ in Eqn. (9) by 
the zeroth order Fy. An approximate form is chosen 

for Sp[^ = 1^0(1911 - kD + - kll) K being 
the total density of the charges and |(7i| a cut-off, corre- 
sponding to the first maximum in Sp{q)). We find : 



ttJ 



In 1 



91 



(10) 



neglecting at first the second term of Sp{q). The 
Coulomb-gas dielectric function of the system is related 
to the charge mobility /ly and to the bound-pair part of 
the dielectric function by |l^] e{uj) — e B + ie^riQ^v iy) / ■ 
We now present the results. Fig. 1 shows, from Eqn. 
(5), 02! (J{)^, as well as (J\I(Jq(^ (that for small a; is 
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essentially fiy (w) the inverse vortex mobility, or vor- 
tex viscosity) versus the logarithmic scaled frequency or 
temperature variable, lny~^. The experimental data 
have been obtained for field-induced free vortices 
(^^^ = 7^ 0) for which regime III is absent 

and g{l) = (1 + {a/^f y^, Va. Thus the coupling 
Ki should scale to zero for / oo (corresponding to 
T > Tkt in the zero field case). We use the simple form 
Ki = KoQilmax - I), and use Imax as fitting parame- 
ter. A good fit is obtained for Imax — 4.81, that is on 
the order of the 1-value for which the linearized scaling 
equations Q yield a vanishing Ki. ai{uj) clearly matches 
Drude behaviour for regime I, 1" 3> 1. The "intermedi- 
ate" regime, Y < 1, with MP dependence '--^ InF is seen 
to be fairly large. The experimental [Q data points for 
SNS arrays, shown in Figs. 1,2, fall in regimes I and II. 
Very low frequency data are not unequivocal and are not 
shown. Typically Q], Fq ^ 300 Hz, uj varies from ^10 Hz 
to - 10 kHz, for Rj ~ 2ml] (SNS arrays) and /./ -100 
nA, and ^(/) ~ 3.1 for / = 0.001, so Y goes from - 0.1 
to - 200. 




0.1 1 10. 100. 



Inverse Frequency Scale (~1/(o^2) 

FIG. 2. Real and, in the inset, imaginary parts of l/e{u}) 
versus an inverse frequency scale ~ y, in a log- log plot. The 
constant-mobility {/ly^luj) = Py) Drude limit is represented 
by the dashed line, while the solid line is the result of our 
calculation b) . The experimental results of Ref. [4] are given 
by the dots. 

Fig. 2 shows the dielectric function e{Lu) (obtained 
through Eqns. (8) and (9), by using the full form of 
Sp{p)) as a function of Y. One again clearly recognizes 
two- frequency regimes (I and II), separated by a crossover 
frequency iOcross ~ ("of^)""^- For lo > LOcross-, 3?(l/e(w)) 
varies like |cli| as in Minnhagen's regime II, whereas for 
UJ < uJcross, it varies like w^, as in Drude's regime I. More- 
over, we have verified for both methods, that the ratio 
Rcr at F = 1 varies between the Drude {Ra = 1) and MP 
{Ra = 2/tt) signatures 

Regime HI is not reached for the F- values shown in Fig. 
1, but for SIS arrays, Rj is orders of magnitude higher, 



so the critical behaviour might be more clearly seen. Be- 
low Tkt, or more generally, for F <C 1, one has effective 
damping coefficients riv — cri due to bound pairs, rather 
than free- vortex inverse mobilities, and CTi ~ l/wZj. This 
is consistent with simulations of driven vortices : there is 
a velocity-dependent viscosity coefficient, decreasing for 
larger velocities Larger oscillating pairs, probed at 
lower UJ, are more sluggish, since they move in a logarith- 
mically interacting viscous medium of smaller pairs. 
In conclusion, we have proposed a regime interpretation 
of anomalous vortex dynamics, based on the ratio of 
the frequency- dependent probe scale, and the screen- 
ing length. Both Drude and anomalous vortex dynamics 
emerge in different regimes, from calculations of the dy- 
namic JJA conductivity and the vortex mobility. This 
reconciles different results, supports postulated conduc- 
tivity scaling, and indicates further dynamical avenues to 
be explored, in simulations and experiments. 
It is a pleasure to thank D. Bormann, P. Martinoli and 
P. Minnhagen for useful conversations, and P. Martinoli 
for reading the manuscript. 
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